Thermal structure of coronal loops as seen with Norikura coronagraph 
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ABSTRACT 

The thermal structure of a coronal loop, both along and across the loop, is vital in 
determining the exact plasma heating mechanism. High resolution spectroscopic obser- 
vations of the off-limb corona were made using 25 cm Norikura coronagraph, located 
at Norikura, Japan. Observations on a number of days, were made simultaneously in 
four forbidden iron emission lines namely, [Fe xi] 7892 A hne, [Fe xiii] 10747 A & 
10798 A lines and [Fe xiv] 5303 A line and on some days made only in [Fe xi] 7892 A 
and [Fe x] 6374 A lines. Using the temperature sensitive emission line ratios [Fe xiv] 
5303 A/[Fe xiii] 10747 A, and [Fe xi] 7892 A/[Fe x] 6374 A we compute the electron 
temperatures along 18 different loop structures observed on different days. We find a 
significant negative temperature gradient in all the structures observed in Fe xiv and 
Fe XIII and a positive temperature gradient in the structures observed in Fe xi and Fe x. 
Combining these results with the previous investigations by Singh and his collaborators, 
we infer that the loop tops, in general, appear hotter when observed in colder lines and 
colder when observed in relatively hotter lines as compared to their coronal foot points. 
We suggest that this contrasting trend observed in the temperature variation along 
the loop structures can be explained by a gradual interaction of different temperature 
plasma. The exact mechanism responsible for this interaction needs further quest and 
has potential to constrain loop heating models. 

Subject headings: Sun: activity — Sun: corona — methods: data analysis — techniques: 
spectroscopic 



1. Introduction 

Plasma in the solar atmosphere is mostly confined to magnetic structures called coronal loops 
which are regarded as the basic building blocks of the solar corona. Any improvement in the knowl- 
edge on these loops, like the exact thermal structure of a coronal loop and other such properties, 



- 2 - 



can therefore help us in solving several long-lasting issues related to corona (for e.g. coronal heat- 



i ng) . One of the earliest attempts to theoretically model the coronal loops, was bv iRosner et al 



(jl978l ). which is widely known as RTV model. They derived scaling laws between temperature, 
pressure and length of the loop and demonstrated that in a u niformly he a ted s table hydrostatic 
loop the temperature maximum must be located near its apex. ISerio et al.l (jl98ll ) generalized this 
model to long loops by including the variation of pressure and heat deposition along the loop. 
Kano &: Tsuneta (jl996l ) derived temperature distribution in 16 steady loops observed with Yohkoh 



in soft X-rays and found that the temperature is highest around the loop top and decreases to- 
wards the foot points consistent with these models. Similar studies on cooler EUV loops, however, 
showed a little or no temperature variation along the loop w hich is not in agreement with the 
static loop models (jLenz et al.l Il999l: lAschwanden et al.l |2000j) . Different theories were proposed 
to explain this isothermal nature. iReale &: PeresI (j200d ) suggested that a superposition of several 
unresolved thi n strands at different teni peratures can produce a flatter temperature profile in the 
observed loop. lAschwanden et al.l (j200ll ) proposed that for long EUV loops mos t of the heating is 



concentrated at the foot points resulting in the near isothermal loop structure. IWinebarger et al 



([2002) showed that inclusion of flows in the coronal loop models (non-static loops) can give similar 
results. Non-uniformity in the loop cross section, with a significant decrease near the foot points can 
also make the temperatu re profile more isoth ermal than that in the case of constant cross section 
(jLandi fc Feldmanlbood ). Some other studies (jPel Zanna Masonlbood : iReale k Ciaravellalbood ) 
disagree on the observed isothermal nature. Later studies mostly focused on the cross-field thermal 
structure of the loops. Here, we investigate the thermal profile along coronal structures using spec- 
troscopic data observed with the 25 cm Norikura coronagraph. We give the observational details 
and data analysis methods in the following sections and discuss the results in the final section. 



Observations 



In view of the complex nature of intensity ratios observed (jSingh et al.ll2004l ). further obser- 
vations were planned and obtained by Jagdev Singh using 25 cm coronagraph for this study. The 
observations were carried out simultaneously in four forbidden iron emission lines namely, [Fe xi] 
7892 A fine, [Fe xiii] 10747 A & 10798 A (IR) hues, and [Fe xiv] 5303 A line, using three different 
CCD cameras (spectra in both the IR lines were imaged on a single CCD). In addition to these, 
some observations were made simultaneously in only tw o lines, [Fe X l ] 789 2 A an d [Fe x] 6374 A. 
The instrumental setup is similar to that described in ISingh et al.l (jl999l . |2003| ) . Four different 
raster scans, with relatively better signal-to-noise, were chosen for this study. Observation time 
and the target region for these four scans are listed in Table [H First three scans were observed 
simultaneously in four lines and the last scan was observed only in two lines. We refer to these 
two sets as set I and set II respectively, for the rest of the paper. The slit length was around 500" 
and the width was 5" for the scans in set I and 3" for those in set II. Each scan covers a portion 
of the equatorial off-limb corona of about 200" x 500" size. Exposure times were varied from day 
to day and also between different lines from 25 s to 50 s depending on the target of observation 
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to get good signal. This translates to a duration of 20 min to 35 min to complete a single scan 
(see Tabled]). Disk spectra were obtained immediately after or before the scan, by keeping the slit 
at the center of the solar disk. Corresponding dark and flat spectra were taken after each set of 
observations. Reference slit images (wire spectra) were also taken with the same setup by keeping 
three wires across the slit at a fixed known distance to correlate the spectral images obtained in 
different lines. 



3. Data Reduction and Analysis 

The data were prepared following the standard reduction procedure for the Norikura datasets. 
Each spectral image had been corrected for the dark current and pixel to pixel variations using 
the corresponding dark and flat images. Scattered light component in the spectra was subtracted 
using the disk spectra. One of the scans in set I (taken on 2005 October 6) was binned on two 
pixels in spectral dimension to improve the signal-to-noise. The dispersion values at each CCD 
was computed by comparing the corresponding disk spectra with the standard solar spectra. The 
dispersion values for the scans in set I are 40.3 mA pixel"^, 74.2 mA pixel"^, and 19.3 mA pixel^^, 
respectively for Fe xi, Fe xiii, and Fe xiv lines (twice the values for the binned case) and the values 
for those in set II are 23.6 mA pixel"^ and 26.6 mA pixel"^ respectively for Fe xi and Fe x lines. All 
the spectra were fitted with a single Gaussian profile and the parameters such as amplitude (peak 
counts), peak position, and line width, were derived. Only those spectra, with a signal-to- noise 5 or 
greater, are considered. Gaussian amplitudes were converted to absolute intensities using the disk 
spectra and standard solar flux values at respective wavelengths. The total area under the fitted 
Gaussian curve is taken as the total intensity in the line profile which is used in this analysis. Peak 
positions were converted to Doppler velocities taking the average position in the scan as reference. 
Li ne widths wer e conv erted to FWHM after correcting for the instrumental broadening as explained 



in lPrasad et al.l ( 20131 ). A 2-D monochromatic image of the observed region was then constructed 



in each of these parameters. Spatial scale of these images in each spectral line is not same due to 



Table 1: Details of observations 



Date 


Target 


Observation 




region 


time (JST^ ) 


2005 Sep 21 


West limb 


09:28 - 09:50 


2005 Oct 3 


East limb 


08:28 - 08:47 


2005 Oct 6 


East limb 


06:05 - 06:41 


2007 Oct 6*^ 


West limb 


11:57 - 12:29 



"Japanese Standard Time 

''Observed only in two lines Fe xi 7892 A and Fe x 6374 A. 
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different focal lengths of the optics used to focus the spectra. This is also affected by the differences 
in pixel sizes. These images were therefore aligned and brought to same spatial scale using the wire 
spectra in each line. The pixel scale after this correction is found to be ~ 2.2" per pixel for the 
scans in set I and ~ 2.4" per pixel for those in set II. However, the spatial resolution in the scan 
direction is limited by the slit width. Final images after all these corrections, for intensities in 
Fe XIV 5303 A, Fe xiii 10747 A & 10798 A, and Fe xi 7892 A emission lines corresponding to the 
scan taken on 2005 September 21, are shown in Figured) 



3.1. Translating Emission Line Intensity Ratios to Temperature 



Ratio of intensities from a suitable line pair can be used to get the thermal information of the 
source. We chose the emission line pairs (Fe xiv 5303 A, Fe xiii 10747 A) and (Fe xi 7892 A, Fe x 
6374 A) to study the thermal structure of a loop, since the emission from both the lines in each 
pair is expected to originate from the same plasma volume due to the closeness in their maximum 
abundance temperatures. The temperatures of maximum abundance for the ions Fe xiv, Fe xiii, 
Fe XI, and Fe x are 2.0 MK, 1.8 MK, 1.4 M K, and 1.1 MK respectively, as noted from the atomic 
database CHIANTI (jK.P. Dere et al.lll997l ). These line ratios are sensitive to electron temperature 
which means that any change in the temperature will be reflected in the ratio of the observed 
intensities. We modele d the dependence of these ratios on temperature using CHIANTI (version 7.1; 
K.P. Dere et"al] (Il997l): iLandi et all koi'i )). Th e calculations are done at a c onstant density using 
the co ronal abundances of ISchmelz et al.l (|2012l ) and ionization fractions of lArnaud &: Raymond 
(|l992l ). Photoexcitation from a radiation field of temperature 6000 K at a distance of 1.1 Rq is 
included. A calibration curve is constructed with the modeled intensity ratio values for a range 
of temperatures. However, it was found that there is a significant dependence of these ratios on 
density as well. We therefore, produced similar calibration curves at different densities varying 
from 10^'^ to lO^''^ cm~^ in steps of 0.1 in the log scale. The const ructed curves for both the ratios 
at four density steps are shown in Figure [2j The calculations by ISrivastava et al.l (|2007l ) indicate 
a weak density dependence of the Fe xiv/Fe xiii line ratio. On the contrary, our analysis using 
CHIANTI shows this to be significant. The density dependence of the Fe Xi/Fe X ratio is relatively 
weaker. Figure [2] also shows an opposite dependence for the two ratios on density. The curves in 
the left panel shift towards the lower temperature side as the density increases whereas that in the 
right panel shift towards the higher temperature side. The temperature at any location of the scan 
can be estimated from these curves, using the density and the observed intensity ratio (from the 
appropriate line pair) at that location. 



3.2. Temperature Variation along Loop Structures 



Individual loop structures that are clearly visible in the images for both the lines of a line 
pair were identified. These structures were followed manually by clicking at many locations along 
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Fe XIV 5303 A Intensity 




100 200 300 400 500 
Distonce olong slit (orcsec) 

Fe XIII 10747 A Intensity 




100 200 300 400 500 
Distonce olong slit (orcsec) 




100 200 300 400 500 
Distonce olong slit (orcsec) 

Fe XI 7892 A Intensity 




100 200 300 400 500 
Distonce olong slit (orcsec) 

Fig. 1. — Intensity images in Fe XIV 5303 A, Fe xiii 10747 A& 10798 A, and Fe xi 7892 A (bottom) 
lines constructed from the scan observed on 2005 September 21. Plus symbols in the top two panels 
locate the individual loop structures chosen. Numbers given at the end of these structures, are for 
reference. 
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the structure in one of the emission hnes and the same pixel locations are automatically chosen 
from the images in other lines observed at the same time. Only those parts of the loops that are 
clearly discernible and seemingly isolated are chosen. Plus symbols in Figure [1] outline the chosen 
loop structures from that scan. Numbers marked at the end of these structures, are for reference. 
Intensity ratios were computed along each of these structures. Electron densities were derived along 
each loop using the observations in IR line pair Fe XIII 10 747 A fc 10798 A whose ratio is sensitive 
to changes in density but a weak function of temperature (jNoens et al.lll984l ). We verified this with 
CHIANTI and found the ratio to be fairly constant with < 1 % changes over a temperature range 
of 1 - 4 MK. However, the signal in 10798 A line is poor and the density information could not be 
obtained at all the locations along a loop. So, we fit the log values of the observed densities linearly 
to derive the value at each pixel location along the structure. Density was found to vary between 
10^'^ to lO^ '^ cm~^ in the chosen structures. Density information could be not obtained for the 
scans in set II as they are observed only in Fe xi and Fe x. For the structures chosen from these 
scans, we model the density variation with the average gradient observed in the structures from 
set I and taking the density at the base to be 10^'^ cm~^. This may have some uncertainty, but at 
least allows us to take the density variation into account. Combining this density information with 
the observed line ratios, we derived the temperature along each loop structure using the calibration 
curves constructed, as explained in § 13.11 These values are then plotted and fitted with a first 
order polynomial to obtain the gradients. We also compute the temperature values for a constant 
density scenario to check the effect of density variation. Observed intensities in 5303 A and 10747 A 
lines, their ratio and the temperature profiles along loop #5 (see Figured]) are shown in Figure [3j 
Temperature values in the bottom left panel are derived for a constant density 10^'^ cm~^ and 
that on the bottom right are derived taking the observed density variation into account. The solid 
lines over plotted on these values represent a linear fit to the data. The scatter in temperature 
values increased when the density dependence of the line ratios is taken into account which could 
possibly arise due to uncertainties in the density estimation. Temperature values at 10" & 100" 
distances from the base of the chosen structure and their gradients are tabulated in Table [2] for all 
the loops investigated in this study. In the last column, we list the temperature gradients in the 
loops keeping the density constant at 10^'^ cm^^, for comparison. The 1-a uncertainties in each of 
these parameters, obtained from the fit, are also given in the table. 



4. Results and Discussion 

We investigated 18 loop structures for the temperature profiles along their length using emission 
line ratios. Out of these 16 are from the scans in set I and the rest two are from those in set II. 
TableOlists the computed temperature at 10" &: 100" distances along with the temperature gradient 
for each of these structures. The gradient values represent a significant change in temperature along 
the loop length (see Figure [3]). Also, none of the structures was observed from the foot point through 
the full length of a loop, so the temperature difference between the foot point and apex of the loop 
can be even larger. Clearly, the gradients are negative for most of the loops (positive but close to 
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Table 2: Temperature values and their gradients along individual loop structures. 



Set 


Date 






loop 


TViTiDPratiirp 


T(^TTi TK^rn tnrp 


V J- d 


yrpb 










(#) 


at 10" (MK) 


at 100" (MK) 


(MK arcsec"^) 


(MK arcsec"^) 


I 


2005 


Sep 


21 


1 


2.12±0.009 


2.05±0.005 


-0.0007±0.00007 


-0.0013±0.00006 




2005 


Sep 


21 


2 


2.23±0.006 


2.02±0.007 


-0.0023±0.00016 


-0.0022±0.00015 




2005 


Sep 


21 


3 


2.26±0.008 


2.04±0.006 


-0.0024±0.00012 


-0.0033±0.00009 




2005 


Sep 


21 


4 


2.02±0.010 


1.92±0.010 


-0.0012±0.00018 


-0.0021±0.00020 




2005 


Sep 


21 


5 


2.05±0.009 


1.97±0.007 


-0.0010±0.00014 


-0.0026±0.00008 




2005 


Sep 


21 


6 


2.05±0.008 


2.07±0.005 


0.0003±0.00011 


-0.0010±0.00009 




2005 


Sep 


21 


7 


2.06±0.008 


2.03±0.005 


-0.0003±0.00010 


-0.0011±0.00007 




2005 


Oct 


3 


8 


2.00±0.014 


1.93±0.009 


-0.0008±0.00018 


-0.0010±0.00014 




2005 


Oct 


3 


9 


2.17±0.013 


2.18±0.008 


0.0001±0.00015 


-0.0007±0.00012 




2005 


Oct 


3 


10 


1.93±0.012 


1.95±0.014 


0.0002±0.00025 


-0.0020±0.00020 




2005 


Oct 


3 


11 


2.30±0.020 


2.29±0.014 


-0.0001±0.00028 


-0.0007±0.00024 




2005 


Oct 


6 


12 


2.17±0.007 


1.95±0.009 


-0.0025±0.00019 


-0.0015±0.00017 




2005 


Oct 


6 


13 


1.99±0.011 


1.97±0.011 


-0.0003±0.00021 


-0.0023±0.00016 




2005 


Oct 


6 


14 


1.99±0.010 


1.79±0.010 


-0.0022±0.00017 


-0.0013±0.00018 




2005 


Oct 


6 


15 


1.96±0.010 


1.92±0.011 


-0.0005±0. 00020 


-0.0014±0.00014 




2005 


Oct 


6 


16 


1.97±0.012 


1.89±0.013 


-0.0009±0. 00022 


-0.0015±0.00016 


II 


2007 Oct 6 


17 


1.21±0.011 


1.49±0.017 


0.0031±0.00022 


0.0034±0.00020 




2007 Oct 6 


18 


1.40±0.012 


1.46±0.017 


0.0006±0.00022 


0.0008±0.00022 



"Temperature gradients taking the density dependence of the hne ratio into account. 
''Temperature gradients for a constant density (10®'^ cm~^). 
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zero for the rest) under set I and positive for those under set II. The average gradient values are 
-0.0009±0.00004 MK arcsec"^ and 0.0018±0. 00016 MK arcsec-^ respectively, for these two cases 
after accounting for the density variation along the loops. Note that the temperature estimation 
in the former case is from [Fe xiv] 5303 A/[Fe xili] 10747 A ratio whereas that in the latter case 
from [Fe xi] 7892 A/[Fe x] 6374 A ratio. The data sample in set II is hmited (only two loops #17 
& #18). However, it complements the extensive analysis done by lSingh et al.l (|2004l ). who finds an 
increase in 7892/6374 line intensity ratio along 90 % of the structures studied, implying a positive 
temperature gradient. Therefore, the observed variations in temperature from [Fe xiv]/[Fe xiii] 
and [Fe xi]/[Fe x], possibly indicate a general behavior. The line pair (Fe xiv, Fe xiii) represents 
a slightly hotter plasma compared to the (Fe xi, Fe x) pair. So, it appears that the loop tops are, 
in general, colder when observed in hotter lines and hotter in relatively colder lines, with respect 
to their coronal foot points. Further, the average temperature values from both the sets are closer 
to each other at larger hei ghts (100") indicating the likelihood to reach a common value at greater 
heights (jSingh et ahlbood ). 



The temperature i n a uniformly heated stable hydrostatic loop increases along its length and 
peaks at the loop top (jRosner et al.lll978l ). Also, if the heating is non- uniform with a bulk of 
the heat deposited close to the foot points, then the temperature maximu m can occur well below 



the loop top leading to a negative temperature gradient along its length (jAntiochos et al. 



199S 



Aschwanden et al.ll200lh . But in such cases, if the heating scale height is sma ller than one-thir d 



of the loo p half-length , loops become thermally unstable and are short lived (jSerio et al.lll98lh . 
Recently, iHuang et al.l ([2012) report the ubiquitous presence of such loops in the low latitude 
quiet corona. However, taking the generality of our results into account, it is not possible to 
explain the contrasting trends observed in different temperature plasma with different types of 
heating. Alternatively, if the loop structures studied here are composed of several unresolved 
strands that are impulsively heated resulting in a multi-thermal structure and if we allow the 
hotter and colder plasma to interact with each other gradually, the observed temperature decrease 
in hotter plasma and increase in colder plasma can be easily explained. Since the loops are observed 
here at a moderate s patial resolution (~ 2" per pixel), it is possible tha t they are multi-stranded 
(jReale Sz Peresll200d ) and multi-thermal (jPatsourakos &: Klimchukll2007l ). At this point, the exact 
mechanism causing the gradual interaction between the hotter and colder plasma remains unclear. 
But we believe that this process holds key to constrain the loop heating models and needs further 
investigation. 



The authors thank the referee for useful comments. J. Singh thanks the staff at Norikura 
Observatory for their help during observations. He also wish to thank NAOJ for providing the 
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George Mason University, the University of Michigan (USA) and the University of Cambridge (UK) . 
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Fig. 2. — Variation of emission line ratios Fe Xiv/Fe Xlll (left) and Fe Xl/Fe X (right) with electron 
temperature modeled using CHIANTI. Different curves are plotted for different densities. See text 
for other parameters used in the model. 
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Fig. 3. — Variation of different parameters along loop #5 marked in Figure [H Temperature values 



in the bottom left plot are derived for a constant density 10 cm~ and that on the bottom right 
are derived taking the density variation into account. The solid lines over plotted on these values 
represent a linear fit. 



